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TRANSFER  FUNCTION,  IMPULSE  RESPONSE  AND 
RERADIATED  WAVEFORM  FOR  AN  ELLIPTICALLY  SYMMETRIC 
RERAD  IATIOII  FUNCTION  OF  THE  FORT-!  **  0 HALF  AN  ODD  INTEGER 
(USL  PROGRAM  NO.  0*37) 

by 

Donald  A.  Stremsky 

USL  Technical  Memorandum  No.  2242-156-67 
1 May  1967 

INTRODUCTION 


A computational  procram  has  been  prepared  by  the  Information 
Processing  Division  to  compute  a particular  Reradiation  Function  v ( X ); 
Transfer  Function  V:  (v,P);  Impulse  Response  w (t,P);  and  Rcradiated 
Waveform  g (t,P)  as  defined  below  in  terms  of  the  incident  plane  'rave 
pulse.  This  IBM  704  programmes  ignated  USL  Program  No.  0837,  is  in 
Fortran  II  language  and  is  described  in  Appendixes  A and  B.  Similar 
computational  programs  are  described  in  USL  Technical  Memorandum  Nos. 
2242-111-67  and  2242-157-67. 

THEORY 

Reference  (a)  contains  a description  of  the  mathematical  model 
constructed  and  the  theory  behind  c inside  ing  reflection  as  a 
reradie.tior.  phenomenon. 

This  program;  computes  for  half-integer  values  of  J) 


This 

program; 

(0 

81  '■ 

(b) 

W (w,P) 

(c) 

Jew  (t,P) 

00 

Kg  (t,P) 

(X) 


j Accrss 
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■■_  ts : program  XA32  as  lescrib  I In  r'3L  T ch.  Mno,  Ho.  224P-157-67 
c.""".  - V'  *.:'■  yc:  ■ . r'.c.i  i'  "i  tirv.s  f-.r  .’nterjer  values  o/V. 


ccitPUTKi-c  rsoc-iui'  D'isc  aFTio;: 

’C'"ic  cliji.  '.,t.r.i)_  f_4-  USI  P ’ : ii  p.  0 37  is  Appenilr.  A, 

• chart  is  / -,.o  ndix  ",  and  the  IIK  7 Oh  Fortran  II  Progjpeu  \s 
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The  basic  input  data  deck  required  by  the  program  consists  of 
four  cards. 


Table  1 
Card  Formats 


Card  No. 


2 


Cols.  Contents 


1-8 

al 

7-16 

a2 

17-24 

X1 

25-32 

x2 

33-40 

c 

4i  -43 

V 

49-51 

V 

52-54 

ISKP  (set  equal  to  zero 

to  compute 

Reradiation  Function) 

55-57 

JSKP  (set  equal  to  zero 

to  compute 

Transfer  Function) 

57-60 

KSKP  (set  equal  to  zero 

to  compute 

Impulse  Response  & reradiated 
waveform) 

For  long  jobs  requiring  the  use  of 
a dump  tape  at  least  one  of  the 
above  variables  should  not  be  set 
equal  to  zero. 

61-63  NSTOP  (in  reference  to  Reradiated 

Waveform  Array  (k,t),  NSTOP  is 
the  numbe  of  times  t is  incremented 
when  k has  its  maximum  value. 


1-8 
9-16 
12-24 
25-32 
35  -36 

39-46 

r 


Initial  value 
Maximum  value 
Initial,  vaoie 
Initial  value 
Kit 


of  U) 

for 
for  t 

of  k (if  not  computed) 


r-i  c 

Av 

7-24 

'5-32 

(if  set  equal  to  zero,  initial 
value  for  k will  be  computed) 

Maximum  value  for  k 

A.  ) 

4 Components  of  f\ 

"?  Components  of  N 
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For rats : 


Cols 

Contents 

33-40 

41-4P. 

49-56 

A Xp 

At 

57-64 

Increment  of  u) 

65-72 

A’* 

1-0 

u)o 

9-16 

£u) 

value  of 

17-24 

Maximum 

25-32 

Initial 

value  of  £ 

33-40 

A© 

4i-48 

f 

6f8.3, 

513 

4f8. 3, 

2X,  12,  2X,  F8.3 

3 - Format  9F8.3 

4 - Format  0FS.3 


Tape  Units  Required 


Tape  Unit  No. 

3 

4 


{ 


5 

6 

7 

r> 

O 

0 

ss 


15  must  be  do  r.  to  dump. 

I.'o  oth.r  sense  switches  are  us«d. 


Tape  Identification 
Data  Input 

Values  for  Reradiation  Function, 
Transfer  Function  & Impulse 
Response. 

Calcomp  Plotter  containing 
va  ues  for  Re radiation  Function 
Reradiated  Waveform  Array  (k,t) 
Transfer  Function  Array  (k,w) 
Impulse  Response  Array  (k,t) 

Dump  Tape 
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Subroutines  Required 

Subroutine  A!.!?  computes  the  va  Lues  for  A array  referred  to  under 
equation  (4). 

Subroutine  SPILT  computes  the  values  of  spherical  Bessel  Functions 
(see  reference  (o)  and  Appendix  C). 

PROGRAM  OUTPUT 

Tape  jjh  contains : 

(1)  The  va  ues  for  the  Z array  plus  the  corresponding  values 
for  the  Reradiation  Function  according  to  Format  (IX,  F1C.5,  5X> 
F10.5). 

(2)  The  values  for  the  product  of  k and  w plus  the  corresponding 
va  ues  of  the  Transfer  Function  according  to  Format  (l)C,  F10.5>  5X} 
F10.5). 

(i)  The  values  for  t/k  plus  the  corresponding  values  for  the 
Impulse  Response  according  to  Format  (IX  F10.5,  SX,  F10.5). 

Tope  ,i 5 contains: 

The  values  for  the  Reradiation  Function  (Calcomp  Plotter  tape). 
Tope  #6  contains: 

The  Reradiateu  Waveform  Array  (k,t)  according  to  Format  (F10.5). 
Tape  $7  contains: 

The  Trai.sfer  Function  Array  (k,w)  w th  Format  (F10.5). 

Tape  contains: 

The  Impulse  Response  Array  according  tc  Format  (F10.5). 
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Tape  .','0  1.  s dump  tape. 

Notes:  This  program  contains  options  to  compute  or  not  to  compute 
a y of  the  functions  mentioned  above.  Tapes  Unit  No's.  6,  7,  and  8 
can  be  used  as  input  to  USL  Program  ,;Od09,  "Representation  of  Surfaces: 
A Computer  Program  to  Plot  Contours  and  Draw  Perspective  Views",  by 
Sdva  d Beardsworth,  Jr. 


SUMMARY 

An  IBM  704  Fortran  program,  USL  Program  No.  0837,  has  been 
written  to  compute  a particular  Reradiation  Function,  Transfer 
Function,  Impulse  Response,  and  Reradiated  Waveform  in  terms  of 
the  incident  plane  wave  pulse. 


<8.  CL. 
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APPENDIX  A 


S (I) 

Z (I) 
RERAD(l) 

TRFER(l) 

AKW(LM.l) 


NOMENCLATURE  LISTING  FOR  USL  PROGRAM  NO.  0837 

Element  of  Reradiation  Function 
Array- 

Element  of  Transfer  Function 
Array 

k \r 


RESP(LM,I) 
RAT  ) 


Element  of  Impulse  Response 
Array 

t/k 


/ 


A1 

A? 

XI 

X2 

C 


Element  of  Reradiated  Waveform 
Array 


x2 


c 


V 

Maximum  value  for  v 
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Initial  value  for  t 


AK 

k 

AKMAX 

Maximum  value  for  k 

B1 

B? 

A,  1 

^ r components  of 

33 

A*»  \ 

B4 

7 components  of 

B5 

A X1 

b6 

A xp 

B? 

At 

Bo 

Increment  of  w 

310 

A^ 

B12 

Ak 

OMEGA 

DELTA 

A 

fjyp 

Maximum  value  of  “V 

TAV 

PHI 

Initial  va  ue  of  "2? 

ib 

appendix  b 

FLOW  CHART  FOR  USL  PROGRAM  NO.  0837 
"RERADIATION  FUNCTION,  TRANSFER  FUNCTION,  IMPULSE  RESPONSE  (CASE  IB)" 
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APPENDIX  C 


FORTRAN  PROGRAM  NO.  0837 


c BfHACtATICN  FLNCTICN. TRANSFER  FUNCT ION  * IMPULSE  RESPONSE  (CASE  lB> 

C C • A # $TPf M SK Y 

DIMENSION  7(500) .RFRAC ( 500) »APPLF(500> *AXW (50*50) «TRFER (50*50) «R6S 
IP  (50*50)  •« AT  IC  ( 50,50)  *0  ( 1 50 ) *S  ( 500 ) *GSUM  (50*50)  ♦«  ( 1000)  * I DUMP  ( 18) 
Cl  PENSION  HuFFfcH(i02A>  *XAX I S ( 500 ) • VAX  IS ( 500 ) 

STCJN  AlF  »0B3? 

STECC  Alf  • 

WRIT*  OUTPUT  TAPE  4*9501 

9501  FCRWAT(IHI) 

BfAC  INPUT  TAPf  3*9502*01 

9502  FcRyAT(A5) 

IF (TCJN-DI >9503*9504 *9503 

9503  PAUSE  ft 

9504  WRlTt  OUTPUT  TAPE  4*9502*01 
WRJTfc  OUTPUT  Tape  4*9505 

9505  FcRPAt(10X32PC.A.STRF^SXY»RCCM  31?6*CO0E  2242) 

READ  INPUT  TAPE  3 . 1 00 . A l * A2 *X l *X? *C »V *N , T5KP * JS<P .KSKP *NSTOP 

100  FCPVAT (6F8. 3*513) 

READ  INPUT  TAPf  3*101  *W«WyAX»T*AX.KIC*A)(MAX 

101  FCRKAT(4F8,3»2X*!2*2X»F8,3) 

Read  INPUT  TAPf  3*103*Bl»e2*B3»e4*B5»B6.P7*B8»B12 

103  FcrvAT(9F8.3> 

READ  INPUT  TAPE  3.104»ANEGA,0ELTA*TT,TAu*B10*PHI 

104  FORMAT (6F8»3> 

Wl»4 

AK1«AK 

ti*t 

TaU1«TAU 

NSTCPsNSTCP*i 

PIE-3.1415 

CEG«180./PIE 

C6«1.?724 

Cl«N*l 

NP1«N*1 

NP2«N*2 

XaN 

C3*X*3*/2* 

C10aX*1.72. 

C?aN*2 

C4«2.««C3 

C5»2.#*C2 

C7*C4*C6/C5 

NPRCCsl 

CC  160  I * l «Np2 
PNEW*?I-1 

PP«CD3NPPCn#NNj;w 

iftc  continue 

PRCDsmPHOD 
IF(ISXP) 116*102*116 

102  1*1 

105  Sm»(Xl/Al>*«2*<X2/A2)«#2 
SX*S(I> 

Z(I)«SORTF(SX> 

ZPlaS(!>-l. 

IF (2P1 ) 110.110.108 
108  RERAD ( I ) *0. 

M*I-1 


CC  TC  261 

11C  *PPLF ( 1 ) * 1 »-S ( I ) 

PEAR*aPPLUP«»C10 

RFPAC(!)sPFA«*C3/PTE 

115  XlsXl*P5 
X2*X2*P6 
I ■ 1 * 1 

GC  TC  105 

261  WPITt  OUTPUT  TAPF  4,251 

251  FORMAT ( 1 X 36HZ  REWA01ATI0N  FUNCTION) 

WRITE  OUTPUT  TAPE  4.252* (ZCI> »«ERAD<T) *I»l*Nl) 

252  FCPN'AT(1X.F10,5.*>X*F10.5) 

WRITE  OUTPUT  TAPF  <,.253 

253  FORMAT!///) 

116  LM*1 

125  IF(JSKP) 195.126.195 

126  Isl 

IF  (KK) 130.12T.130 

127  Pl»fll-C»B3/V 
P2«e2-C«6<*/V 

RBs(At«Pl>##2* (A2«P2)«»2 
AK*SCPTF(RB>/C 
130  Isl 

194  WRITE  OUTPUT  TAPE  4.254 

254  FORMAT (1X34HKW  TRANSFEP  FUNCTION) 

165  IF(AK-AKVAX) 170,170*195 

170  IF(W-WMAX)  180, 160.190 
160  AKW<LM»I)sAK*W 
IF  (W)  165*182*185 
182  TRFER(LM»I)*1.00n00 
GC  TC  187 
185  BCPsAKW(LM.I) 

CALL  SPNJ(BOP»B) 

BANG»PCP«*C3 

TRFER(LM.l )sC7»PRCD«B(NP2)/BANG 
187  1=1*1 
WsW*B8 
GC  TC  170 

190  N2*I-1 

N4sLM 

WRlTt  'OUTPUT  TAPF  4.255* ( (AKW(LM,I>  »TRFER <LM • I ) * la  1 »N2) ,LM«N4,N4) 

255  FCRMAT(1X*F10,5*5X.F10.5) 

280  bRITt  OUTPUT  TAPE  7 ,273 » ( (TRFER (LM* I ) , Is l »N2) *LN«N4.N4) 

273  FORMAT (F10. 5) 

1*1 
W*W  1 
LM=LM*1 
AK* A<  *F  1 2 
GC  TC  165 

195  1*1 

end  file  7 
END  FILE  7 
LM*1 

ak*aki 

1F(KSKP)250.265,250 
265  t»R  I Tfc  OUTPUT  TAPF  4,257 

257  format «ix34mt/k  impulse  response) 

191  ApSFTrABSF (T> 


n 


IF  (APSFT-AK)  199*199 #193 
193  RESP  <L^* I » 

N'3*I-1 
GO  TC  235 

199  NPRCCal 

200  CC  220  J»1»NP1 
NE*J*J 

NPRC0«NPRCD#NEa'J 

220  continue 

PRCC2«NPRCR 
RATICaM»I)«T/AK 
62»PATI0(LW«I)*»2 
G3«1,.G2 
G4*G3**C1 
CEN«C5*PRCC2 
RESP  (LN  ♦ I ) *PRCD*'%*  /DEN 
T»T*R7 
1*1*1 
GC  TC  191 

235  N3«I-J 
N4*Lm 

kRlTt  OUTPUT  TA^E  4 , 258  » , (RATIO (LM •I)*RESP(tM«I)*I*l*N3) #LM*N4*N4) 
2*8  FCRPATtlX.Fi0.bt«i*:.FlC.5) 

282  WRITE  OUTPUT  TAPE  8 .2 75 ♦ ( <RE5P «LM* I) . I *1 *N3 ) *LM*N4 ,N4) 

275  FORMAT (F10.5> 

N3Pl«N3*l 

IF(NSTCP-N3Pl)23A,23W»23l 

231  CC  232  I«N3Pl  .t:sT0P 
RESP (NA.l) *0*0 

232  CONTINUE 

WRITE  OUTPUT  TAPF  8,233*  < <RESP<LM«n  «UN3Pi»N$T0P>  *LM»N4*N4> 

233  format  ifio.5) 

234  !F(AK-AKMAX)236'237«237 

236  Ak*AK*P12 
tM«LM*l 
1*1 

T*T1 

GO  TC  199 

237  Lm*  l 

EnC  FJLE  8 
ENC  FILE  8 
N5*0 
Ak»AX1 
1*1 

238  Cll«B7/AK 
T»T  1 

J*  1 

Cl2*Afj5F  (Oil) 

239  G5UV(LP»J)«0* 

240  FRACT«TAL7TT 
El*AMpGA*CELTA«FWACT/2#0 
E2*E l*T  AL 

E3«E?*PHI 
E4«CrsF (E3/DEG) 

CALL  ANP(TAU*R) 

FCN*R»!)*E4 
TCIF* (T-TAUJ /A* 

CRAPE*ABSF (TCIF) 


1 6 


PLUN*GPAPE-1 .0 
IF  ( P L t r.‘ > 8 5 o • e 5 r.  • 2 r «> 

279  E9*0.P 

C-C  TC  F01 
950  E5*GPAPE/CL2 
NE5«E5 
lA«NE5*l 
IP*NE3*2 

Efc*ORAPE-HAT!C(L^»!A> 

E7*f 6/T  ll 
Efl«l.r-E7 

BSPNSf  7*«ESP<L“.IP) ♦E8#RE5P  <l*«IAt 
E9«FCv*wsPNi#en 
901  C5lN<Lf'«J>*GSL  '(I.M.JWE9 
TflL*TAl ♦PlO 
I«T*l 

E10«AP5F ( T AU  > 

IF(TT.EU)>241  .240.249 

241  G5GMLN'*J)»GSL^(LF'.J>/AK 
AbSFT  * ABSF ( T ) 

IF  (AHSM-AK)242,243.243 

242  T,T*P7 
J*J*1 
Tau«tAU1 

IF ( SENSE  SwITCm  5)9997.9999 

9997  CC  9999  LtC«l«15 
ICUNP(UO**0 

9998  CONTINUE 
ICUPP ( 16) *-6 
IDUNP  < 1 7 ) =*0 
ICUMPUB)  *N5 
CALL  DL'MP(IDUNP) 

9999  GC  TC  239 

243  IF  |AK-AKNAX)244.244»2A9 

244  AK«AK*B12 
N3«J-1 
N4«LV 

276  WRITE  OUTPUT  TaPF  6 .277. ( (GSUM (LM»1 ) , I *) »N3) »LM»N4  »N4) 

277  FCBNAT (F10.5) 

N3Pl*\3*l 

IF (KSTCP-N3P1) ^24 .224.221 

221  DC  222  I«N3Pl.NSTnP 
GSLP(NA.I) *0.00000 

222  CONTINUE 

WRlTt  OUTPUT  TAPE  6 .223 • ( (GSUM <LM « I) « I «N3P l .N5TCP ) »LM«N4 ,N4 ) 

223  FCRVAj  ( F 10 . 5 ) 

224  N5*N5^J^NSTOp-N3 

IF (AK-AKNAX) 2 IR, 2 78 .249 

278  LV*L^*: 

T *T  1 

TaUsTaUI 

1*1 

C-C  TC  2 38 

249  E*C  FILE  6 

end  file  6 

250  N i *M ♦ i 
N2«Ni*2 

283  IF (I5<P) ?50.259,150 


I*) 


259  CALL  PL0TS(BUFFER(102A) *102A«5) 

DC  2 60  J=l.Nl 

XAXISIJ )*Z(J) 

YAXlS(J) sRERAD  < J) 

260  CONTINUE 

CALL  PLOT  (0*0.5.0*-3) 

CALL  SCALE  <YAXIS»5.0«Nl*l»10.0) 

CALL  SCALE  <XAXI5«10.0»M*1*10.C) 

CALL  LINE  (XAX!S*YAXIS.N1«1*1*11) 

CALL  AXIS  <0*0.0.0.20hRE«ADIATICN  FUNCTICN«20*5.0f90.0*YAXI5<Mi > *V 

lAXIStw?) *10.0) 

CALL  AXIS  (0.0.0.0»1HZ»-1.10.0*C.P*XAXI5(^1) »XAXIS(W2) *10.0) 

CALL  PLOT  (0*0.0.0*999) 

350  READ  INPUT  TAPE  3, 9502. ED 
IF  (ED-TECD)9503. 9509*9503 

9509  WRlTt  OUTPUT  TAPE  A. 9511 
9511  FCRPAT(AHCENO) 

END  FILE  A 

9510  STOP  5 

END (1*1 *0*1*1) 
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c 

c 


SPHERICAL  BESSEL  FUNCTION  J 


R.  0.  WHITTAKER 


0614 


SUBROUTINE  SPHJ  (x.B) 

DIMENSION  B (200) »BL (200) «SC (4) »A (2) 
BE*l.OF-10 
CC  5 1=1.200 
PL  < I > =0.0 

5 P(I)=0.0 
XsXM.C 

S STC  A (2) 

S STC  XL 

S STC  A (1) 

IF  (X-.05)  45.6.6 

6 IF  (X-100.)  7.7.45 

7 IF  (X-10.)  0.15.15 

8 RN«72./(4.02-L0GF(X)I 
C-C  TO  20 

15  RN=1.51*X*25. 

20  HAXsRN*  1.0 
NsHAX-2 

B (MAX-l ) =BE*1 .0 
S STC  Bl(MAX-l) 

CC  35  1 = 1. N 
JsMAX* I ♦ 1 
AN* ( 2* J-3 ) 

AN*AN* 1.0 
CALL  DPA1 
S CIA  B(J-l) 

5 LOG  BL(J-l) 

CALL  DPA13 
S STC  VH 

S STC  VL 

5 CLA  X 

S LOG  XL 

CALL  PPA1 
S CLA  VH 

S LOG  VL 

CALL  DPA  14 
S STO  VH 

S STC  VL 

S CL*  B('j) 

S LOG  BL(J) 

CALL  DPA 1 
S CLA  VH 

S LOG  VL 

CALL  DPA  1 2 
S STC  B ( J-2 ) 

5 STC  HLU-2) 

35  CONTINUE 

CALL  CPSC  (A (2) »SC (4) . ICUHHY) 

S CLA  X 

S LOG  XL 

CALL  DPA1 
S CLA  P(l) 

S LDC  Bl(l) 

CALL  DPA  l 3 
CALL  DPA 1 
5 CLA  SC (4 ) 


2.1 


